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Abstract 


A  method  to  determine  the  variance  of  statistical  fluctuations  for  instnanents 
used  in  the  registration  of  the  cosmic  radiation  -  the  neutron  monitor,  the  duplex 
cubical  counter  telescope  and  the  directional  telescopes  -  is  presented.  The 
standard  errors  of  data  from  cosmic  ray  instruments  are  estimated.  The  results 
are  compared  with  the  standard  errors  calculated  from  the  theoretical  Poisson 
distribution.  It  is  stressed  that  the  Poisson  distribution  vill  give  an  under¬ 
estimation  of  the  standrad  errors. 
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System  of  notations 

The  following  system  of  notations  is  used: 

In  the  population: 

D(x)  The  standard  deviation  of  the  variable  x 
2 

D  (x)  The  variance  of  the  variable  x 

E(x)  The  mean  of  the  variable  x 

C(xy)  The  covariance  of  the  variables  x  and  y. 

In  the  sample: 

2 

S  The  variance  calculated  from  the  sample 

S(x)  The  standard  error  of  the  variable  x 

X  The  mean  of  the  variable  x 

^  indicates  estimation 
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I.  Introduction 


In  dealing  with  investigationt:  where  counting  tubes  are  used,  it 
to  assume  that  each  observation  is  dravm  from  a  Poisson  distribution, 
known  that  a  series  of  random  events  will  give  the  probability  that  n 
occur  in  the  interval  t  (for  instance  Cun’an  &  Graggs,  19^9)* 

/  ^  .  \n  -  Art 

P(x  .  i 


is  common 
It  is  well- 
events  will 


(1) 


where  A  is  the  average  nxomber  of  events  per  second.  As  long  as  the  events  will 
occxir  at  random,  which  holds  for  most  radioactive  processes  and  for  cosmic 
radiation,  the  Poisson  distribution  will  be  theoretically  exact.  This  problem  has 
been  theoretically  discussed  by  several  authors  (Bateman  1910  and  Fry  1928).  In 
connection  with  Bateman's  paper  of  I9IO,  Rutherford  &  Geiger  (I9IO)  made  their 
famous  experiment  on  the  distribution  of  o(  -particles  from  a  radioactive  sample. 
The  result  showed  good  agi'eement  with  the  Poisson  distribution. 

This  gives  the  variance  of  the  counting  rate  N 

(X)  =  N  (2) 

The  dlsciate  Poisson  distribution  called  "the  law  of  small  numbers",  is  skew 
for  small  values  of  At  but  when  Increasing  Xt  the  distribution  forms  more 
symmetrical  and  will  soon  be  approximately  equal  to  the  continuous  normal 
distribution.  According  to  Wolfenden  (19^1)  the  ske\mes8  of  the  distribution  will 
be  negligible  for 

At  =  N  =  10 


In  registrations  of  cosmic  radiation  the  use  of  counting  tubes  are  common  and 
for  error  calculations,  the  Poisson  distribution  is  of  basic  Importance.  One 
must  be  aware  that  this  estimate  will  only  take  care  of  the  true  random  variations 
in  the  number  of  pai-ticles.  It  is  sometimes  used  uncritically,  regardless  of  the 
existence  of  other  types  of  errors  which  may  be  introduced  by  the  Instzamients  and 
during  the  data  processing. 

In  this  paper  the  standard  errors  of  cosmic  ray  data  from  different  instruments 
are  cedculated  and  the  results  are  discussed  and  compared  with  expected  theoretical 
values . 


II.  Statistical  treatment 


There  are  difficulties  in  error  calculations  of  time  series  aneGyses.  The 
genezal  formula  for  the estimation  of  the  variance  is: 


D^(x) 


1 

n-1 


n 


If  the  variable  x  is  a  function  of  time,  the  result  will  be  an  estimate 
including  both  the  time  variations  of  the  variable  and  the  statistical  random 
fluctuations.  Thus  this  method  will  be  unsuitable  for  our  purpose. 

Me  Cracken  (1958)  has  chosen  a  very  calm  cosmic  radiation  period  for  such 
a  oedculatlon  and  got  a  fairly  good  estimate  .-f  the  variance  T^{x)  due  to  only 
statistical  fluctuations,  although  the  existence  of  a  daily  variation  and  a 
possible  27-day  variation  will  result  in  overestimation  of  ri2(x). 

To  pass  these  problems  there  opens  a  possibility  by  the  IGY  recommendations 
to  divide  the  cosmic  ray  instruments  in  equal  sections  to  get  an  unbroken 
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continuous  registration.  By  serving  each  identical  section  of  the  instrument 
through  a  complete  set  of  electronics  a  good  defence  against  an  all-over  break¬ 
down  is  offered.  Each  section  is  working  to  some  extent  as  an  independent 
instrument.  Only  the  power  line  is  common.  The  sections  are  measuring  the  same 
phenomena.  By  calculating  the  differences  or  the  ratifs  between  those  two  sections 
we  have  a  method  to  get  the  statistical  random  fluctuations. 

Each  observed  value  can  be  written  as 


N  =  E(I)  +  f^,(P,T)  +  fg  +  f  (5) 

=  the  true  mean  of  the  intensity  of  the  cosmic  radiation  during  the 
period  of  observation 

=  the  atmospheric  influence  of  the  cosmic  radiation  'r-he3E  the  pressure 
P  and  the  temperature  T  are  functions  of  time> 

=  the  time  variation  of  the  primary  radiation.  We  assume  that  the 
variation  at  the  top  of  the  atmosphere  is  proportional  to  the 
variation  at  the  sea  level  for  constant  atmospheric  parameters 
(Ehmert). 

£.  =  random  variable  due  to  statistical  fluctuations,  which  will  consist 
of  both  the  Poisson  error  and  the  instrumental  errors. 

We  know  that  E  (  t*  )  =  0 


where  E(l) 
f^  (P,T) 
fg  (  isl) 


There  exist  different  theories  of  the  atmospheric  influence  of  the  radiation 
(Dorman,  Duperier  etc.)  but  for  this  calculation  it  is  svifficient  to  use  the 
simple  linear  form  of  the  wellknown  Duperier 's  formula 

dN  =:  <<  ap  +  /S'dT  +  xdH  (4) 

dN  is  the  difference  between  the  observation  and  the  corrected  vsdue;  dP,  dT  and 
dH  are  the  differences  between  the  means  and  the  observed  values  of  pressure, 
temperature  and  height  respectively,  to  a  certain  pressure  level  (100  or  200  mb)j 
oi  ,  and  'S  are  constants. 

Our  knowledge  of  the  time  variations  of  the  primary  radiation  is  limited.  We 
know  to  some  extent  the  periodical  variations  such  as  those  with  periods  of  2k 
hours  and  27  days  but  the  nonperiodicial  ones,  i.e.  Forbush  decreases,  eure  still 
uncertain. 

The  two  identically  built  sections  record  the  same  radiation.  We  Introduce 
the  ratio  between  the  true  counting  rates: 


(5) 


where  I  and  I  are  the  number  of  counts  in  the  i*esp.  sections  due  to  cosmic 

*  y 

radiation  when  the  atmosphere  is  supposed  to  be  constant  from  time  to  time. 

z  is  the  number  of  covmts  recorded  in  one  section  that  due  to  multiplicity  in 
the  Instrument  will  be  detected  in  the  adjacent-  sectlbn  simultaneously. 

The  time  variables  P,  T,  H,  z  and  I  ai’e  not  true  random.  For  a  longer  period 
the  distributions  of  the  variables  tend  ?o  be  symetrical  around  the  mean.  P,  T  and 
I  are  independent  and  we  assume  that  their  covariances  will  be  negligible.  This 


will  not  be  too  misinterpreted,  especially  as  we  in  the  follwoing  aro  using  the 
differences  of  (x-y),  which  will  reduce  the  influence  of  existing  covariance 
term.  But 

z  .  f^ 


AH  =  f^  (AP) 


-k. 


In  the  following  the  symbols  D,  E  and  C  will  mesui  the  varianAe,  mean  €uid 
covariance  indei)endent  of  the  variables  are  true  random  or  not. 

We  know 

E(Ci.P)  =  E(/iT)  =  E(  AH)  =  0 

The  registrated  values  from  each  section  are  divided  by  using  eg.  (J),  (^*-)  and 
(5)  for  a  certain  time 


*1  “ 

I  +  z,  +(I  + 
^  1 

=  kl 

yi 

+  z  +(kl  +  z. 
-L  y^ 

*«: 

(-■ 

11 

I  +  Zt  +  (I 

^  ^1 

we  have: 

E(x) 

=  kE(l  )  +  E(z) 
«y 

E(y) 

=  E(I  )  +  E(z) 

«/ 

We  form: 

A  oO‘‘  ^D^(P)  +  /'^ 

(6) 

(7) 


and 


D^(x)»k^D^(l  )+D^(z)+a|  kV(l  )+E^(2)+kV(l  )+3^(z)j  +2k  C  (l  z)(i+a| 


1 


+  2kD^(P)E(l  )  E(z)  +  D^(  .1  ) 

/  A 


D^(y)=D^(lJ+I'^(z)+A  [e^CI  )+E^(z)+D^(l  )+D^{z)l  +2C(l  z)  (  1+a) 
y  ^y  j  y*!^^ 

^  A 


+  2  D^(P)  E(y  E(z)  +  D^(  £y) 


c(xy)  =  kD^(l  )  +D^(a)  +A  i  kE^(l  )+E^(z)+kD^(l  )+D^(z)]  + 
[]c(lyZ)  Jl+Aj  +D^(P)E  (ly)  E(z)]  (k+1) 

Difference:  The  differences  (x-y)  give  from  eq.  (6)  and  (7): 


(8) 


(9) 


(10) 


E(x-y)  =  (k-1)  E(I  )  =  M 
D^(x-y)  =  D^(x)  +  D^(y)  -  2  C(xy) 

where 

C(xy)  /  0 

as  X  md  y  are  correlated.  Eq.  (8),  (9)  and  (10)  give: 
D^(x-y)=(k-l)^ 

’*'”,,,,2  (Ed,)  -  2(1^))  ^  ^ 

S^(I  )  E^( I  ) 

'  y  y 


(11) 


-5- 


We  are  interested  in  the  distributions  of  C  ^  and  £  ,  which  describe  the  true 
random  fluctuations.  *  ^ 


Ratio:  We  form 

X,  .  kl  ,  +  +  (kl  A  Pj  gA%)+  f, 


“  1  +  +n  +  z^ O ■  P“ + /5 « +  yajr^')+ f y 


The  variance  of  a  function  F(x^  x^  -  -  -  )  can  be  formed  (KendaJJ.) 

1:1^=  <*1^)} 

when  i  /  k 

{^)  =  -5^  2  - -^g'^  ■■  [E(xy)-E(x)  E(y)3 

y  E^Cy)  E^y)  E(y)E2(y)  ^ 


If  I  and  kl  >>  z  and  ki=s  1  we  can  write: 

y  y 


0  [  D^(z)+  a(e^(z)+  D^(z)}1  +  [d^(£  )+kV(  e  )\  (12) 

y 


E"(ly) 


2 

If  we  assume  z  =  D  (z)  =  0 

(5)  «  ■  I  (  £  )  +  (  £  )1 

'y  E^(i)  L.  '  x'  '  y'J 

^  of  Uj)^ 

Rema£ks:Due  to  eq.  (11)  the  variance ^^a~dependent  of  k,  the  atmospheric  effects 
Md  the  variations  of  the  primary  radiation. 

The  distribution  ol  the  ratio  can  be  assumed  approximately  normal  under  the 
following  conditions: 


(13) 


E(I  )  and  E(I  )  >>  0 

X  y 

I  and  I  £  and  £  respectively. 

X  y  X  y 

This  is  fulfilled  in  most  cases  when  we  are  dealing  with  observations  of 
cosmic  radiation.  The  variance  is  due  to  eq.  (12)  Independent  of  atmospheric  and 
eosmical  variations  for  z  =  0. 


III.  Discussion  of  errors  introdi^ed  by  the  method  of  registration 

In  registration  of  cosmic  radiation  it  is  often  useful  to  scale  the  data 
before  the  recording.  The  registered  value  is  sometimes  rounded  off  to  make  the 
following  data  processing  as  simple  as  possible.  The  factor  used  in  scaling  and 
rounding  off  is  chosen  according  to  the  counting  rate.  We  shall  here  discuss  the 
errors  introduced  by  the  methods. 

Scal^i^:  Scaling  can  be  made  in  different  ways  but  at  continuous  registration  it 
is  a  rule  to  scale  without  zero  resetting.  This  means  that  all  counts  ai’rlved  in 
the  scaler  and  not  scaled  is  carried  over  to  next  period  of  registration. 

We  assume  at  the  time  of  registration  a  fraction  of  m,  the  scaling  factor, 
consisting  of  N  counts  has  arrived  in  the  scaler  and  will  not  be  registered.  The 
following  registration  will  then  start  with  N  counts.  For  one  registration  zre 
write  ^ 

X,  =  mN,  -i-  -  N 

1  1  f  e 

where  mN^^  is  the  number  of  counts,  which  are  registered  as  scaled.  N  and  can  be 
assumed  ‘‘"as  independent.  We  form  ® 

E(xJ  «  nE(Nj  +  E(N,)  -  S(N  )  (iK^ 

life  \  / 


-6- 


D^(x^)  =  n^(N^)  +  (Ng  +  N^)  +  2  mC  [n^,  (N^  +  *\)J  (15) 

To  draw  conclusions  of  the  corrections,  which  must  be  applied  to  the  variance  of 
the  scaled  distribution  to  get  the  true  distribution,  is  rather  complicated.  She 
influence  of  the  covariance  term  is  depeiident  of  the  shape  of  the  distribution, 
the  magnitude  of  the  sca].ing  factor  and  how  the  group  intei’vals  ai’e  situated  in 
the  distribution.  As  long  as  the  scaling  unit  is  chosen  small  compai'ed  with  the 
standard  deviation  of  the  distribution  of  x  the  Infliience  can  be  negligib3.e.  In 
this  paper  all  observations  aice  scaled  with  a  small  unit  and  no  corrections  ax*e 
made. 

The  distributions  of  N  and  N  are  discrete  with  the  range  (0,m).  If  the 
distribution  of  x  is  appx’oxfmately^symetrical,,  which  is  the  case  when  the  period 
of  observation  is  lorg,  we  can  assume  that 

E(Ng)  =  E(N^). 

Thus  the  method  will  not  introduce  any  systematical  error  according  to  eq-  (l4), 

Roundlng_cff_^  When  tabulating  data  it  is  often  useable  to  round  off  the  data  to 
simplify  the  foUwoing  data  process.  Eq.  (1^)  will  then  be  formed  when  q  is  the 
rounding  off  factor 

E(x^)  =  qE(N^)  +  E(N^)  (l6) 

The  method  of  roxonding  off  is  chosen  that  E(N  )  is  as  close  to  zero  as  possib.le. 
One  method  of  roui\dir.g  off  often  used  is,  thaf  all  values  ended  with 


are  increased- with  one  vunit.  E(N  )  is  Aot  exactly  zero  but  the  diffeience  is  small,. 

e 

The  distribution  of.  N  is  discre-te  xrtth . the .  range 

e 

(-  f»  +  f) 

In  this  paper  no  rounded  off  values  ax:e  used. 

Influence  of  the  recording  time.  As  yet  xre  have  assumed  the  time  of  the  recording 
as  sero.  It  is  desii-able  to  hold  this  time  as  short  as  posslb3.e.  Photographic 
methods  are  then  excellent  hut  mechanical  printing  counters  are  difflcu3.t  to  bxiild 
with  short  printing  time.  Durii-xg  the  time  rf  recording  process  therc'-ia-a  protthl'llty 
that  ■.  counts  has  arrived  in  the  scal.er.  The  mechanical  counter  is  not  open  for 
counts  during  printing  and  a  unit  is  in  this  case  missed.  This  wi3.1  i’esult,  that 
E(N  )  <  E(N  )  in  eq.  (l4)  and  a  systematic  error  is  introduced.  As  a  rule  the 
tln«  of  priSting  is  small  e.g.  the  mechanical  counters  used  for  the  counter 
telescopes  at  Uppsala  and  Murchison  Bay  have  a  printing  time  of  0.8  sec.  For  jthe 
international  cube,  which  has  a  counting  rate  of  each  section  of  about  5  x  10 
counts/hour,  every  tenth  recording  ■wi.ll  ralSs  one  unit.' The  introduced -ryst'^natica]. 
error -Wil'l  be  of  tv.-s  magntuUde  of  5  per  cent  of  m. 

W  Resultsj_ 

Neutron  monitor:  In  order  to  get  cn  estimation  o'i  the  standared  deviation  of 
the  true  statisticeul.  f luctuat.  i  ons  for  an  IGY  standard  neutron  monitor  (Gimpson 
1953)  data  from  three  different  stations  have  been  xised. 


lat. 
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long. 

Altitude 

Mean  counting  rate 

Murchison  Bay 

80°03'N 

18°15  'E 

Sea  level 

per  hour 

26000 

Uppsala 

59°51'N 

17°55'E 

Sea  level 

25000 

Mt.  Wellington 

42°55 'S 

147°i4  'E 

725  m 

18000 

The  stations  are  all  equipped  with  neutron  monitors  built  after  the  IGY 
recommendations  in  two  sections,  each  served  by  a  complete  set  of  electronics. 
Great  care  has  been  taken  to  use  such  periods  for  the  calculation,  'vrtien  the  two 
sections  of  the  Instruments  have  been  working  without  cheinge  of  connting.  Diagram 
1  shows  5?  for  10-day  periods  from  the  three  stations  during  the  periods  used  in 
the  calculations. 

The  differences  between  the  two  sections  are  calculated  for  one-hour  values 
uncorrected  for  atmospheric  effects.  We  then  have  the  sample  mean  and  variance  by 
the  usual  statistical  formulae: 

d  =  i  /I  (x-y) 
n 

S^(x-y)  =  i  I(d-d)^ 

where  d  is  the  difference.  We  estimate 
d  •=  E  (x-y) 

S^(x-y)  =  ^  (x-y) 

For  large  values  of  n  we  find  that  is  a  good  estimate  of  D^.  This  happens 
in  most  cases  in  this  paper. 


The  standard  errors  are  calculated  from 

S(d)  = 

(17) 

Vn 

S(s2)  =  ^\"2n 

(18) 

S(S)  . 

(19) 

These  equations  are  valid  when  the  variable  is 

drawn  from  a  normal 

distribution.  Diagram  (2)  shows  the  distribution  of  (x-y)  for  Murchison  Bay  with 
the  normal  distribution  fitted.  Ab  seen  the  distribution  can  be  dasumed 
approximately  normal.  „ 

Eq.  (11)  shows  that  calculating  D‘^(x-y)  the  result  will  be  an  overestimate 
of  p  p 

D^(£^)  +  D^(ey)  -  ' 

for  M  0.  We  can  calculate  the  degree  of  overestimation  if  we  know  M,  E(I  ), 
D^Iy),  d2(P),  d2(t),  d2(H),o<,/I,  and  TJ  .  ^ 

For  the  neutron  monitor: 

o  o 

D  (P)  and  D  (I  )  are  more  difficult  to  estimate.  Meteorological  records  give 
information  of^D2(p).  For  Stockholm,  which  is  a  typical  sea  level  station  at  a 
middle-high  latitude,  a  mean  covering  10  years  gives 

"D  (P)  -  12  mb 

This  value  can  be  expected  to  Increase  during  simmer  and  decrease  during  winter. 
The  variations  are  small. 
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O 

D  (I  )  is  the  variance  of  the  primary  radiation  and  vill  vary  greatly 
according^to  whether  the  i)eriod  choosen  for  the  calculation  is  cosmically  calm  or 
not  and  at  which  latitude  the  station  is  situated.  Diagram  3  is  drawn  for  the 
stations  at  Uppsala  and  Murchison  Bay  for  different  vedues  of  d2  (I  )  and  D'-(P). 
Following  assumptions  are  made:  ^ 

E(l  )  >=  12500  counts/hour 

y 

=  0.73  i/mb 

D(I  )  =  100,  500,  1000  counts/hour 

y 

D(P)  =  9,  12,  15  mb 

From  the  diagram  3  It  is  easily  shown  that  for  small  values  of  M  the 
overestimation  will  be  negligible.  If  we  can  accept  an  overestimate  of  1  per  cent 

of  o  2 

D  (  £  )  +  (  e  ) 

X  y 

(for  Uppsala  and  Murchison  Bay  400  counts/hour)  we  can  use  M  —  200  counts/hour. 

Table  number  1  shows  D^(x-y)  for  separate  neutron  monitors.  The  result  shows 
poor  agreement  with  the  theoretical  value  given  by  the  Poisson  distribution. 
According  to  Me  Cracken  (1958)  this  is  due  to  multiplicity  of  the  particle 
production  in  the  atmosphere  as  well  as  to  star  production  in  the  moiiitor  itself. 

He  has  shown  theoretically  how  the  multiplicity  gives  a  higher  value  of  the 
standard  deviation  than  expected  from  the  Poisson  distribution.  He  used  experimental 
veilues  from  a  small  two-counter  monitor  and  foimd  the  standard  deviation  of  the 
distribution:  _ 

B(n)  =  1.13  Vn 

irtiere  N  is  the  counting  rate.  He  made  a  rough  estimate  for  a  standnrd  I.G.Y. 
monitor  conslting  of  12  counters 
D(N)  ^1.2  VI 

A  mean  covering  three  periods  from  Murchison  Bay  and  Uppsala  which  have  monitors 
of  similar  appearance  and  electronics  give 
D(N)  =  1,181  '/N 

The  result  shows  good  agreement  with  the  one  of  Me  Cracken.  The  result  from  Mt. 

Wellington  can  not  be  Immediately  compared  with  this  vedue  as  *=  46?  counts/hour. 

We  can  then  use  the  ratio  between  one  hour  observations  of  the  two  sections 

to  find  the  distribution  of  6  and  £  .  By  calciLlating  «)  we  get  an  over- 

X  y  j 

estimate  of  _ 

from  eq.  (12)  due  to  z.  low  we  know  that  z  is  only  a  small  fraction  of  the  counting 
rate  according  to  Me  Cracken  (I958), 

We  assume: 

z  ~  4  per  cent  of  the  counting  rate 
D(P)  10  mb 

'^O.J  per  cent  per  mb 
D(z)  ) 

\dilch  give  a^negllgible  overestimate.  Eq.  (I3)  can  be  vused. 

Table  1  gives  the  result  from  Mt.  Wellington  for  the  calculation  of  the 

ratio 

D(N)  »  (1.25  *  0.02)  V'N 
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Thls  value  is  markedly  higher  than  those  from  Uppsala  and  Murchison  Bay.  The  data 
of  Mt.  Wellington  are  scaled  with  a  factor  6k  while  the  Uppsala  and  Murchison  Bay 
data  are  unsealed.  The  correction  of  the  "scaled"  variance  to  get  the  "true" 
variance  is  difficult  to  find  as  seen  in  a  preceding  section.  However^  the 
correction  is  small.  Shepeard's  correc-^ion  for  the  variance  -l/l2  (see  Cramer 
19**-5)  gives  a  rough  estimate. 

The  difference  between  the  calculated  variances  indicates  that  D(N)  is  an 
instrumental  constant.  The  method  used  to  estimate  the  variance,  will  include  «.n 
statistical  errors  introduced  by  the  method  of  registration,  the  Poisson  error  as 
well  as  the  Instrumental  errors.  The  neutron  monitors  may  differ  in  appearance 
and  electronics,  which  can  cause  a  variability  of  D(N)  from  one  station  to 
another  and  of  cotirse  also  from  time  to  time  according  to  the  condition  of  the 
instrument. 

Counter  telescope.  The  variance  of  the  statistical  fluctuations  for  the  duplex 
cubical  counter  telescope  and  the  directional  telescope  are  calculated  by  the  same 
method  as  in  the  section  of  the  neutron  monitor.  2-hour  values  uncorrected  for 
atmospheric  effects  from  the  two  stations  at  Uppsala  and  Murchison  Bay,  which  are 
equipped  with  a  duplex  cubical  counter  telescope  built  after  I.G.Y.  recommendations 
and  a  telescope  for  east  and  west  directions,  are  used.  Descriptions  of  the 
instruments  are  found  elsewhere  (Sandstrbm,  Dyring,  Lindgren  i960).  As  in  the 
calculation  of  the  neutron  monitor  great  care  is  here  taken  to  chose  periods  of 
observations  where  no  changes  in  the  counting  rate  occur.  It  is  sometimes 
difficult  to  find  such  reliable  periods  for  counter  telescopes  due  to  their 
construction. 

Duplex  cubical  counter  telescope.  Fig.  shows  the  daily  mean  intensity  during 
the  periods  used  in  the  calculation  for  the  international  cube.  If  there  exists  a 
difference  between  the  counting  rates  of  the  two  sections  D2(x-y)  there  is  an 
overestimation  of  the  sum  of  (  £  )  and  (  £  )  acc.  to  eq.  (11).  For  the 
duplex  cubical  counter  telescope  we*aseume:  ^ 

=  0,15  per  cent/mb 

mb  level) 
mb  level) 


mb  level) 
mb  level) 

Correlation  coeff.  (PH)  «  0.5 
E(ly,)=  10^  coiints/  2  houi’ 

D(I  )=  5000  counte/2  hour 

These  values  are  roughly  estimated.  For  a  maximal  overeat imat ion  of  1  per  cent 
we  can  allow 

M  —900  counts/  2  hour 

The  result  of  the  calculations  are  shown  in  table  2.  The  mean  over  all  periods  for 
each  station  gives  ^  _ 

Uppdala  p(N)  =  I.II8  VN 

Murchison  Bay  D(N)  =  1.277  ^ 

Directional  telesc^ie.The  Instruments  at  Uppsala  and  Murchison  Bay  have  eight 
channels,  four  in  each  direction.  The  channels  directed  east  have  odd  and  the 
west  one  even  numbers.  Figs.  5,  6,  and  7  shew  the  daily  mean  intensity  during  the 
periods  used  in  the  calculations.  We  make  roughly  the  same  assumptions  as  for  the 
international  cube  a  part  from 

D(ly^)  =  3*10  counts/2  hour 

D(N)  =  1500  counts/  2  hour 


/^=  0.1  per  cent/°C 

(200 

y  =  6  per  cent/km 

(100 

D(P)  =  10  mb 

D(T)  =  10°C 

(200 

D(H)  =0.5  kmv 

(100 

10- 


Then  we  may  allow  for  2 -hour  values 
M  -s  500  coimts 


for  a  maximum  overestimation  of  1  per  cent.  The  four  channels  of  each  direction 
give  six  sets  of  differences.  OnD.y  two  are  Independent.  Thus  we  have  six  equations 
and  four  variables  £  >  C  )  E  ,  iL  -  hy  the  common  least  square  method  the  variables 

X  y  2  V 

are  estimated.  The  results  are  shown  in  tables  3  11''?  mean  over  all  periods 

for  each  station  give 


Uppsala  p(K)  =  1.102  /F 
Murchison  E0yD(N)  =  I.I67  /n 


Remarh£.  The  influence  of  the  multiplicity  (Me  Cracken  1958)  is  smaller  for  a  counter 
telescope  than  for  a  neutron  monitor.  The  time  dealy  between  mesons  produced  by  the 
same  primary  pai'ticle,  which  can  be  detected  in  each  counter  telescope  section, 
will  mostly  be  small.  The  resolving  time  is  not  sufficient  to  separate  particles 
derived  from  the  same  primary.  The  local  production  of  mesons  is  also  small.  Thus 
we  can  expect  a  closer  agreement  to  the  Poisson  distribution  for  a  counter  telescope 
than  for  a  neutron  monitor.  Still  the  results  do  not  Indicate  this  and  we  find 
further  a  difference  in  the  variance  of  from  Uppsala  and  Murchison  Bay  and  between 
different  channels  of  the  same  instrument  as  well  as  from  period  to  period.  This 
indicates  that  the  standard  deviation  of  £  is  an  instrumental  constant.  The  counter 
telescope,  which  is  equipped  with  GM-counters,  univibratoi's,  shar-peners,  coincidence 
circuits, and  scalers.  Introduces  errors  due  to  resolving  time,  spurious  counts  etc., 
and  are  comparatively  sensitive  to  changes  in  the  power  line  and  trimming  conditions. 
These  types  of  errors  are  more  prominent  on  a  counter  teD.escops  than  on  the  neutron 
monitor. 

The  station  of  Murchison  Bay  was  situated  in  an  area  with  hax-d  weather 
conditions  and  its  electric  source  was  diesel-engined  generators.  Compared  to  the 
more  quiet  conditions  of  Uppsala  it  is  not  sui'prising  to  find  a  higher  variance  of 
€  for  Murchison  Bay. 

It  must  be  stressed  that  the  periods  of  calculations  are  very  carefully 
chosen.  As  a  rule  one  is  forced  to  accept  data  with  less  accuracy  to  get  contlnous 
registrations.  Thus  it  is  obvious  that  the  Poisson  distribution  a  counter  telescopes 
gives  an  undei'estimation  in  calculating  the  standard  errors  of  the  data. 

Pressure  of  the  atmosphere.  To  correct  the  cosmic  ray  data  for  atmospheric  effects 
it  is  necessa::y  to  have  continuous  registration  of  the  atmospheric  pressure.  At  the 
cosmic  ray  station  of  Uppsala  this  is  made  by  a  precision  aneroid  barometer,  which 
is  photographed  simultaneously  ^rith  the  cosmic  ray  recording.  The  mean  pressure  of 
2-hour  periods  is  calculated  from 

P  ri(P„  .  P^  +  ?,J  i 

where  P^,  P^  and  P^  are  the  measurements  at  the  start,  in  the  middle  and  at  the  end 
of  the  2-ho\ir  period  respectively. 

At  Murchison  Bay  a  Plf-hour  barograph  was  used,  which  every  third  hour  was 
calibrated  by  a  standard  mercurial  barometer.  The  2-ho\ir  mean  was  taken  directly 
from  the  barograph  registration-  This  method  is  better  than  the  miethod  used  at  Upp¬ 
sala  ^en  the  pressure  variations  ere  non-linear.  However  as  a  rule  the  difference 
between  the  methods  are  very  small. 

To  get  an  estimate  of  the  variance  of  £  for  pressure  data  we  have  cedculated 
the  difference  in  recordings  fi’om  the  Uppsala  Cosmic  Ray  Station  and  the  Meteorolo¬ 
gical  Institution  of  Uppsala  University.  The  two  statiens  are  lying  about  5  km  apa-.*t 
smd  at  the  same  height  above  the  sea  D.evel,  We  assume  that  the  stations  registrate 
the  same  pressuvre.  The  data  from  the  Meteorological  Institution  is  determined  from  a 
24-hour  barogi'aph  which  is  corrected  by  readings  of  a  good  mercurial  barometer  three 
times  a  day  at  03.10,  14.10,  and  20.10  L.T..  In  the  calculation  cnl;/  the  veLLues  at 
08, l4  and  20  L.T,  have  been  used  to  get  the  best  determined  obsei'’.'aticns .  The 
registrations  are  independent  and  we  have 
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f  c>  *  ”"  <  *«> 

where  C  indicates  the  cosmic  ray  station  and  M  the  meteorological  station.  Fig.  8 
shows  the  monthly  mean  and  the  variance  of  'the  year  Oct.  1959  “  Sept. 

i960.  For  the  whole  year  we  have 

D(  +  D(  =  0.300  -  0.007  mb 

1  The  method  of  recording  the  atmospheri-  pressure  is  supposed  to  be  more  accurate 

j  at  the  meteorological  station.  We  then  estimate 

j  D(f^)  =  0.2mb 

As  a  rule  we  want  the  mean  pressure  over  a  period  e.g.  2  hours.  The  standard 
error  of  such  a  pressure  data  will  be  dependent  of  how  many  observations  are  made, 
i  the  method  used  and  the  time  variations  of  the  pressure  during  the  period. 

( 

I  Conclusions 

1  Data_uncorTected_for_atmos£her^c_effects_^  Our  calculations  have  given  estimates  of 

the  error  variances  from  different  cosmic  ray  instruments.  As  long  as  we  choose 
periods  during  which  the  instruments  are  working  reliably,  £  is  a  random  variable. 

'  D( £  )  is  then  the  standard  deviation  of  the  £  -distribution  at  a  given  tjrue  counting 

rate  N,  equal  the  mean  intensity  of  cosmic  radiation  during  the  period  of  observation 
By  taking  the  ratio  between  the  calculated  standard  deviation  and  the  standard 
f  deviation  estimated  from  the  Poisson  distribution  we  got  D( £  )  as  a  function  of  the 

!  counting  rate  N.  For  simplicity  we  write  the  moments  of  the  t  -distribution  in 

this  section  as  E(n),  D(T),  D^(P)  and  so  on,  since  the  parameters  N,  T,  P...  are 
assumed  to  have  only  random  statistical  fluctuations. 

The  standard  deviations  in  this  paper  are  calculated  from  periods  during  which 
the  instruments  have  been  working  satisfactorily.  From  our  results  an  estimate  of  the 
magnitude  of  D(N)  can  be  made  for  uncorrected  data  from  Uppsala  and  Murchison  Bay 

IGY  neutron  monitor:  D(N)  =  1.2  Vn 
Counter  telescopes:  6^(N)  =  I.I5VN 

It  has  already  been  pointed  out  that  these  figures  might  differ  from  time  to 
'  time  between  different  instruments. 

'  Data  corrected  for  atmospheric  effect^.  The  registrated  cosmic  ray  intensity  at 

ground  levels  Ts  a  function  of  atmospheric  parameters-  To  use  these  data  in  time 
series  analyses  they  have  to  be  corrected.  In  this  paper  we  use  simple  linear 
Duperler  formula  and  at  a  certain  time  we  have 

where  N  is  the  corrected  and  N  the  registrated  value;  P  ,  T  and  H..  are  the 
observeS  values  and  P  ,  T  and  5  are  mean  values.  In  the ‘^following  T’,  T  and  H 
indicates  the  differences'^etween  the  observed  values  and  the  means.  The  error 
varaiance  of  is 

D^(N^)=D^(N^)  +  D^(N^P  cC  )+D^(N^T^)+D^(N^H'>f)+2D^(N^)  B  (20) 

where  ->  .  _ 

B  =E(r)E(P)  [i+E(/9)E(T)  +  E(  ■iJ)E(H)J  +  E(p)E(T)  |_1+E(  «  )  E{H)J  +E(  )E(H) 

and  the  variance  of  a  product  abc  is 

D^(abc)==E(bc)D^(a)  +E(ac)D®(b)+r(ab)D^(c)+D^(a)D^(b)D^(c) 

’.'•i  have  a  .sumed  all  variables  as  independent.  This  is  not  exactly  true  for  P  and  H 

buL,  their  covariance  term  may  be  neglected.  In  the  first  approximation  we  may  also 
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2  2  2 

neglect  D  (<<),  D(^),  D(?j)  and  their  products.  We  further  use  the  following 
error  variances  and  constants: 


D^(N  )  =  1.11  N 
'  r  r 

D^(N  )  *  1,30  t 
r  r 

D®{P)  =  O.Ol  mh^ 
D^(T)  =  5(°C)^ 
‘D^(H)  =  25x10"^  km^ 


IGY  neutron  monitor 
Counter  telescopes 


(from  this  paper) 

(  ) 

(from  Trefall,  NordS,  1959) 
(  ) 

(  ) 


^  •-0.73  per  cent/mb  IGY  neutron  monitor 

if 

P  =  +0.1  per  cent/  C  Counter  telescopes 
1<  =  -6  per  cent/km  -"- 

For  long  periods  when  we  can  assume 
P  =  T  =  H  =  0 

the  error  variance  of  corrected  value  is 

/vO 

IGY  neutron  monitor:  D(N)=l.llK  + 

c  r 

yvg  . 

Counter  telescopes:  D  (N^)  «  1.30  + 

As  seen  from  eq.(2l)  this  is  not  valid  when  P,  T  and  H  differ  from  zero. 

D(N  )  is  dependent  of  the  true  observed  value  of  the  atmospheric  parameters.  This 
is  ^important  when  we  want  to  calculate  the  standard  errors  of  short  period  data 
e.g.  2-hour  values.  Eq.  (21)  will  then  give  the  error  variance.  However,  eq.(20) 
and  (22)  will  in  many  cases  give  a  good  estimate.  The  time  variations  of  the 
atmospheric  parameters  will  affect  the  standard  errors  of  the  corrected  data 
according  to  during  which  period  the  means  of  the  parameters  are  calculated. 


2x10’^  (21) 

r 

ll.lxlO"^  (22) 


Reme^k£:  It  must  be  pointed  out  the  i'.iTportance  of  estimating  the  standard  errors  of 
data  with  care  due  to  which  kind  of  analysis  they  will  be  used  for.  For  fine 
structure  interpreting  it  might  sometimes  be  necessary  to  make  carefxil  calculations 
but  on  the  other  side  rough  estimates  are  often  satisfying.  The  labour  of  the 
error  calculations  must  be  weighted  against  the  need  of  accurate  estimations. 
However,  it  must  be  stressed  that  the  use  of  the  Poisson  distribution  gives  an 
underestimation  of  the  standard  errors  of  cosmic  ray  data. 
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Errata 

In  Fig.  8  is  written  S(Pj,-  P^^).  Read  D(Pj^-Pj,^). 
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Table  ^4-  a 


Dire£WoMl  telescope  a't_l^P£ala_ 


12/8  -  l/ll  1957 
5/11  -  9/12  1957 
25/12-  9/2  1957-58 
12/3  -  '^/*»  1958 
10/6  -  27/9  1958 


f  21/8  -  15/10  1957 

g  5/11  -  9/12  1957 

h  25/12-  25/1  1957-58 

i  51/7  -  1/9  1958 


Channel 

Period 

£  ) 

acc.  to 
eq.  (2) 

S 

S 

P 

§(|) 

P 

1 

a 

32288 

28747 

1.16 

0.03 

b 

36560 

28226 

1.14 

0.05  : 

c 

33628 

29221 

1.07 

0.05  1 

d 

32329 

28729 

1.06 

0.05 

e 

35858 

28426 

1.12 

0.03 

2 

f 

517I+I 

28603 

1.05 

0.04 

g 

41529 

28205 

1.21 

0.06 

h 

55975 

29354 

1.58 

0.07  ' 

i 

35063 

29012 

1.10 

0.05  I 

3 

a 

31649 

29206 

1.04 

0.05 

b 

35970 

28835 

1.12 

0.05 

c 

53716 

29805 

1.34 

0.06 

d 

56751 

29389 

1.39 

0.07 

e 

27737 

29000 

1.00 

0.03 

f 

32391 

28234 

1  07 

0.04 

g 

23360 

27846 

0.93 

0.05 

h 

41147 

28833 

1.20 

0.06 

■ 

i 

28568 

28577 

1.00 

0.05 

5 

a 

39959 

28411 

1.19 

0.04 

b 

30186 

28153 

1.05 

0.04 

c 

31317 

29226 

1.03 

0.04  1 

( 

6 

f 

33302 

28407 

1.08 

0.04 

g 

32774 

28122 

1.08 

0.05 

h 

27679 

29540 

1.08 

0.05 

i 

52614 

28255 

1.08 

0.05 

7 

a 

41574 

29022  1 

1.20 

0.04 

b 

27629 

28714  1 

0.98 

0.05 

c 

29236 

29784  1 

0.99 

0.04 

■ 

d 

50082 

28735  ! 

1.02 

0.05 

e 

55277 

28447  I 

1.11 

0.03 

8 

f 

57582 

29185  j 

1.12 

0.04 

g 

55408 

28605  1 

1.11 

0.05  ; 

, 

h 

36558 

29904  1 

1.10  j 

0.05  I 

i 

-j 

41592 

28217  j 

1.21 

0.06 

Table  ^4-  b 


Directional  telescope  at  Murchison_Bay_ 


a  i4/9  -  lO/lO  1957  d  15/9  -  10/10  1957 

b  29/12-  21/2  1957-58  e  4/1  -  28/2  1958 

c  12/5  -  1V6  1958  f  18/4  -  17/5  1958 


Channel 

Period 

=  D^( e  ) 

p  P 
acc.  to 
eq.  (2) 

S 

P 

t(|  ) 

P 

1 

a 

55301 

29336 

1.10 

0.06 

b 

31752 

30501 

1.02 

0.04 

c 

32004 

29211 

1.05 

0.05 

2 

d 

35811 

29622 

1.10 

0.C6 

e 

55285 

30452 

1.35 

0.05 

f 

40575 

89404 

1.18 

0.06 

5 

a 

5059^ 

28750 

1.33 

0.07 

b 

3*4683 

29797 

1.08 

0.04 

c 

41812 

28711 

1.21 

0.06 

4 

d 

55182 

28942 

1.38 

0.08 

e 

32515 

29840 

1.04 

0.04 

5 

a 

36558 

29055 

1.12 

0.06 

b 

40654 

29955 

1.16 

0.05 

c 

47283 

28856 

1.28 

0.06 

6 

d 

44431 

28716 

1.24 

0.07 

e 

36245 

29581 

1.11 

0.04 

f 

*45925 

28485 

1.27 

0.07 

7 

a 

28470 

29388 

0.99 

0.05 

b 

56978 

30198 

1.37 

0.05 

8 

d 

47228 

28870 

1.28 

0.07 

e 

49124 

29751 

1.29 

0.05 

f 

36011 

28772 

1.12 

0.06 

10-day  means  of  (x—y)  from  the  neutroBi-jaonitors* 


Number  of  events 


•M  H  Vt  'H 


Fig. k.  Dally  mean  intaoBlty  for  the  duhlex  cubical  counter  taleacopes 
st  U^eala  and  Nurchlaor.  Bay. 


Coit'dirMtiwi  of  U’trojf  Mtteop*  Uppoola 


MM  H  MM 

Pig. 7*  Daily  wan  intensity  for  the  directional  telescope  at  Uppsala. 
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